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ABSTRACT 

The atmospheric water vapor content above the Roque de los Muchachos Obser- 
vatory (ORM) obtained from Global Positioning Systems (GPS) is presented. GPS 
measurements have been evaluated by comparison with 940nm-radiometer observa- 
tions. Statistical analysis of GPS measurements points to ORM as an observing site 
with suitable conditions for infrared (IR) observations, with a median column of pre- 
cipitable water vapor (PWV) of 3.8 mm. PWV presents a clear seasonal behavior, 
being Winter and Spring the best seasons for IR observations. The percentage of 
nighttime showing PWV values smaller than 3 mm is over 60% in February, March 
and April. We have also estimated the temporal variability of water vapor content 
at the ORM. A summary of PWV statistical results at different astronomical sites is 
presented, recalling that these values are not directly comparable as a result of the 
differences in the techniques used to recorded the data. 

Key words: Site Testing — Infrared: General — Instrumentation: miscellaneous 



1 INTRODUCTION 

The triatomic molecules H2O, CO2 and O3 are the main 
responsible for reducing the few transparent windows on 
the infrared atmospheric transmission spectrum, and pro- 
ducing absorption bands that are difficult to correct during 
the processing of astronomical data. Out of these windows, 
the atmosphere is opaque and observations at those wave- 
lengths have to be carried out from space. Sites with large 
atmospheric water vapor levels have narrow and unstable 
transparent windows. High levels of water vapor reduce the 
atmospheric transparency, but also increase the thermal in- 
frared background. 

It is assumed among the astronomical community that 
sites placed at lower altitudes are less suitable for high- 
quality IR observations; however, other parameters such 
the tropos phere thickness might be al so playing an impor- 
tant role (|Garci'a-Lorenzo et al.|[2004 ). Indeed, the evalu- 
ation of the IR quality at the Observatorio del Teide in 
the Canary Islands, Spain (| Mountain et al. 1985h and the 
comparison with prediction models ( Cohen et al.lll993 ) in- 
dicate that El Teide site, at ~2400 m above sea level, 
might be as good as Mauna Kea (4200 m) for IR astron- 
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omy in the 1-5 /im window (|Hammerslevl 1 19981 ). Obser- 
vational results pointing to a similar conclusion were re- 
ported for Roque de los Muchachos observatory (ORM) 
on the island of La Palma (Canary Islands, Spain), 
also at ^--^2400 m above sea level, when comparing in- 
frared observations from the 3.58m Telecopio Nazionale 
Galileo with data obtained at the 10m Keck telescope 
on Mauna Kea astronomical observatory (Hawaii, USA) 



(http://www.tng.iac.es/news/2003/03/21/nics_refurbishl 



There are many parameters accounting for the quality 
of an astronomical site, namely seeing, cloud cover, ground 
winds, high-altitude winds, etc (see iMufioz-Tufion et al.l 
(2007) for a review on the characterization of these parame- 
ters at ORM) . The water vapor content is an important pa- 
rameter affecting the IR quality of astronomical sites. Con- 
ditions for IR astronomical observations were classified in 
terms of precipitable water vapor (PWV hereafter) in four 
divisions (|Kidger et al.ll 19981 ): 1) Good or excellent — > PWV 
^ 3 mm; 2) Fair or mediocre -> 3 < PWV ^ 6 mm; 3) 
Poor — >• 6 < PWV < 10 mm; and (4) Extremely poor — >• 
PWV ^ 10 mm . The development of IR instrumentation 
and the requirements for current large and future extremely 
large telescopes demand a proper characterization of PWV 
and statistical studies of large temporal databases (covering 
years). The fraction of nights with good IR conditions (small 
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column of water vapor) as a function of the epoch of the 
year will allow an optimal scheduling of telescope observing 
time. The total atmospheric water vapor contained in a ver- 
tical column of unit cross-sectional area extending between 
any two specified levels is known as precipitable water va- 
por and it is commonly expressed in terms of the height to 
which that water substance would stand if completely con- 
dense d and collected in a vessel of the same unit cross sec- 
tion (| American Meteorological Society! |2000| ). PWV is also 
referr ed to as the total column water vapor (|Ferrare et alj 
|2002| ). Measurements of PWV can be obtained in a number 
of ways, from in situ measurements (radiosondes) to remote 
sensing techniques (photometers, radiometers, GPS, Imag- 
ing Spectroradiometers on satellites, etc). Radiosondes have 
been the primary tn sttu observing system for monitoring 
water vapor; however, the operability of radiosondes is lim- 
ited due to running cos ts and decreasi ng sensor performance 
in cold dry conditions (|Li et al.ll2003l ). Usually, radiosondes 
are expected to measure PWV with an uncertainty of a few 
millimetres, which is considered to be the stand ard accu- 
racy of PWV for meteorologists HieiretaLlllool). Near in- 
frared (NIR) radiometers measuring methodology assumes 
plane-parallel atmosphere, hence it is only satisfactory for 
measurements near the Zenith. Moreover, their uses are lim- 
ited to photometric and bright nights (from first quarter to 
last quarter Moon). Errors of 20% or below are estimated 
for PWV measurements using NIR radiometers in the range 
3-10 /xm, whe reas errors in th e 0.5-1/im range might be as 
high as 40% (|Quesa dal ll989l ). Microwave water vapor ra- 
diometers (WVR) observations are not reliable when liquid 
wate r is present on the WVR frequency window (|Liou et al.l 
I2OOII ) or when there is significant scattering fro m liquid wa- 
ter d roplets and ice crystals in the field-of-view (jZhang et al.l 
Il999l ). 

GPS is an increasingly useful tool for measuring PWV, 
which has gained a lot of attention in the meteorological 
community. The GPS procedure to estimate the PWV is 
based on the fact that the propagation of electromagnetic 
waves through the atmosphere is drastically affected by vari- 
ations on the refraction index of the troposphere, which de- 
pends on the water vapor pressure, air pressure, and tem- 
perature. The consequence is an induced delay in any signal 
crossing the atmosphere, being due to a combination of a 
hydrostatic and a water vapor delay. The hydrostatic delay 
is very stable and has a direct relationship with local at- 
mospheric pressure. The second component, the wet delay, 
is directly related to the water vapor content in the atmo- 
sphere above the site where the measurements are taken. 
The rapid temporal variations of the water vapor affects 
its prediction and proper measurement. The overall tropo- 
spheric delay at a GPS station allows the estimation of the 
PWV with a high-deg ree of accuran cy (iBevis et al.lll994 
iBoccolari et a l. 2006; Jin et al.ll2009l l. although under very 
dry e nvironments GPS mig ht underestimate the PWV con- 
tent (jSchneider et al.ll2009f ). GPS system provides a better 
spatial coverage and continuo us PWV estim ations in com- 
parison with other techniques l|Ge et al.''2000'') . The compar- 
ison of GPS, radiosondes and W VR (Emardso n et al.ll2000l : 
iNiell et al.ll200ll : iLi eral] 120031 ) shows agreement generally 
at the level of 1-2 mm of PWV, corresponding to 7-13 mm 
of zenith wet delay (error << 10%). 

In this paper, we present statistical results on the PWV 



above the Roque de los Muchachos Observatory (La Palma, 
Spain) derived from GPS measurements for the period span- 
ning from June 2001 to December 2008. For comparison pur- 
poses, we also present PWV statistical results for Mauna 
Kea site derived also from GPS data. 



2 DATA SOURCES 

2.1 PWV measurements at ORM 

The Roque de los Muchachos Observatory is located at an 
altitude of 2396 m above sea level on the island of La Palma 
(Canary Islands, Spain) at latitude 28° 45 North and longi- 
tude 17" 53 West. This observatory is one of the final candi- 
date sites for the location of the 42 m European Extremely 
Large Telescope. The Spanish Instituto Geografico Nacional 
has installed permanently since May 2001 a GPS receiver at 
ORM. The GPS antenna is placed on a very stable monu- 
ment which ensures that the GPS satellite data recorded at 
the station does not contain any spurious antenna movement 
that could affect the scientific exploitation of the GPS data. 
This permanent GPS is part of the international network - 
European Reference Organisation for Quality Assured Breast 
Screenin g and Diagnostic Services (ww w.euref.eu) - of GPS 
stations (|Kenveres fc Bruvninxl [20041 ) . The data from the 
GPS station at the ORM for the period spanning from June 
2001 to December 2008 were processed by Soluciones Avan- 
zadas Canarias S.L. (www.sacsl.es), a company devoted to 
GPS data processing for t he International Glob al Navigation 
Satellite Systems Service l|Romero et al.ll2003l ). Global GPS 
data were processed using the dual frequency data recorded 
at around 35 stations worldwide including the station at 
ORM. The daily data processing runs take data from each 
station using all the data available above 10 degrees of ele- 
vation at each ground station every 5 minutes (epoch). The 
satellite orbits are considered fixed as downloaded from the 
International Global Navigation Satellite Systems Service 
(www.igs.org). The Tropospheric Zenith Delay (TZD) for 
each station is estimated every two hours with an initial 
value of 2.1 meters and an a-priori sigma of 2 meters, hence 
being considered "free" . The long standing stations have co- 
ordinates/velocities p ublished in the Intern ational Terres- 
trial Reference Frame (lAltamimi et al.l l2007') which are very 
precise and can be kept almost fixed. The coordinates for 
the stations in study are estimated using the initial condi- 
tions with a 3 mm a-priori sigma. This allows the coordinate 
estimation to be adjusted slightly over the processing arc to 
accommodate small daily variations from the published In- 
ternational terrestrial reference frame positions. 

In each observation epoch, each station has measure- 
ments to all the satellites in-view. Nominally each station 
records measurements to 8-12 satellites per epoch. This 
means that for the TZD estimation (one value every two 
hours per station) each TZD value is averaged over 24 epochs 
(5-minutes each). Therefore, more than ^ 200 individual 
measurements (24 epochs x 8-12 satellites) contribute to 
the TZD estimation. A two-hour averaging time is standard 
in GPS processing for TZD, although shorter estimations 
can be performed. More details on the theoretical basis for 
the ground-based GPS technique, the backgro und physics 
and potential applications can be found in e.g. iBevis et al.l 
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(| 19921 ). The accuracy of GPS measurements are not limited 
by the system receiver noise, but largely by modeling er- 
rors. The uncertainties of the system are multidimensional, 
including errors from satellite orbits, errors in the mapping 
functions used for the atmosphere delay estimations, scat- 
tering in the vicinity of the GPS antenna, the effect of any 
radome in the antenna, etc. An introduction and descrip- 
tion of GPS d ata proce s sing c an be found, for example, in 
iBlewitd l| 19981 ) or lLeickl l|2004l ). 

The TZD derived from the GPS data processing is com- 
posed of two contributions, the hydrostatic (ZHD) and the 
wet (ZWD) delays. The ZHD can be calculated through 
the latitude (4> in degrees), altitude (H in meters) of the 
GPS station (|Saastamoinenlll973l . ll972l ') and the atmospheric 
pressure (P in hPa) according to the following expression: 



ZHD 



22.768 X 10"" X P 



26.6 X lO-4cos(20) - 0.28 x lO-^ x H 



(1) 



The ZWD accounts for the PWV in the atmosphere 
and can be derived by subtracting ZHD from the TZD 
value. PWV can be then calculated as a proportion 
of the estimated ZWD through the following equation 
l|Askne fc Nordiuslll987l ): 



PWV = ZWD X 



10" 



pRv[{k3/Tm) + k2] 



(2) 



where p is the density of liquid water, Kv is the specific 
gas constant for water vapor, k2=70.4 K mbar~ ^, k3=3.739 
lO'^ K^mbar-i are the refra ctivity constants JSevis et al.l 
1 1994 lAskne fc Nordiusll 19871 ) , which have been determined 
by direct measur ements using microwave techniques (see 
iBevis et all (|l994l ) for more details on these constants). Tm 
is a weighted mean temperature of the atmosphere, which is 
strongly correlated wi th surface tempera ture (Ta) according 
to the linear relation (|Bevis et al.lll993 ): 



70.2 + 0.72T, 



(3) 



The uncertainties in the refractivity coefficients of the 
water vapor and the selected mapping functions in the TZD 
calculation can introduce a scale factor in t he PWV derived 
from GPS measurements (|Niell et al.ll200lf ). Therefore, the 
comparison of PWV estimations from GPS data with other 
techniques will require an initial model consisting of a scale 
factor between techniques and an offset accounting for a 
bias in the measurements. This bias comes from residual 
delays that are estimated as corrections to a delay model 
assuming a standard atmosphere using a zenith angle de- 
pendent mapping function (Klciier 2001). The estimation 
of the PWV from TZD derived by GPS measurements re- 
quires surface pressure and temperature measurements. The 
ZWD will not include an additional bias if accurate surface 
pressure is used in the estimation of ZHD. Unfortunately, 
there is not a weather station at the pr ecise locatio n o f the 
GPS antenna at ORM; hence a model l|Boehm et ar.ll2007l ) 
of the pressure and temperature at the longitude, latitude 
and height of the station for each date was used to obtain 
PWV from the derived TZD using equation [T] [2 andO The 
GPS time series for ORM consist of more than 31500 indi- 
vidual estimations with a two-hour temporal resolution of 
the PWV for the period June 2001 to December 2008. 

During the site prospection for the emplacement of the 
10 meters Gran Telescopic de Canarias, different campaigns 



of PWV measurements above ORM were carried out from 
February 2000 to July 2002 using a two band radiometer. 
The radiomet er operates in the 940-nm wa t er v apor band 
(|Kidger et al.| [i998. 200j; IPiniUa et al1l2002l : IPinilla _2003i). 
Observations were taken on the Moon (solar observations 
are also possible including an appropiated neutral density 
filter) , measuring the entire 940-nm water vapor absorption 
line (filter centered on 946.7 nm with a bandwidth of 19.8 
nm) and nearby continuum (filter centered on 884.5 nm with 
a bandwidth of 18.2 nm). The PWV is directly correlated 
to the square of the ratio of the signal in the two filters. 
No model-dependent assumptions are made for this calcula- 
tion. The absolute calibration of the 940nm-radiometer was 
measured from radiosonde flights for a range spanning from 
0.5 to 25 mm of water vapor. This remote sensing tech- 
nique assumes a plane-parallel atmosphere, restricting the 
measurements t o nearly Zenithal positions ((Ouesada 1981; 
iGao et al.l[l993 ) in order to convert the measured line-of- 
sight PWV to the zenithal values. The relative precision of 
the line-of-sight PW V is around 5% for a single measurement 
(jKideer et al.lll998l ). The two main drawbacks of this tech- 
nique are that is only usable in photometric conditions and 
with the Moon well above the horizon. 940nm-Radiometer 
required a human supervision and the temporal resolution 
was about 30 minutes. 

A direct comparison between GPS and the 940nm- 
radiometer could be carried out for the period from June 
2001 to July 2002. In this sense, 243 PWV estimations 
from the 940nm-radiometer were recorded. In order to per- 
form the comparison, we re-sampled the 940nm-radiometer 
dataset to the same temporal resolution of the GPS data. 
Following the temporal re-sampling, the number of data was 
reduced to 55 in both the GPS and the radiometer. As we 
already mentioned, errors up to 40 % has been reported for 
PWV estimations w hen using radio meter measurements in 
the 0.5-1 fim range (jOuesada Il989l ). We have adopted the 
standard deviation of the averaged 940nm-radiometer data 
in the two-hour sampling as the error of the PWV estima- 
tion, being smaller than 25% for all the estimations. Errors 
of about 1 mm are estimated for the PWV derived from the 
proccesing of the GPS data. We obtained a good correla- 
tion between GPS and 940nm-radiometer measurements of 
PWV for the period June 2001- July 2002 (Fig. [ij, with a 
Pearson correlation coefficient of 0.97. A linear relationship 
was determined using the task LINMIX_ERR ( ,Kellv 2007i ). 
available from the IDL Astronomy User's Library, which in- 
corporates a Bayesian approach to linear regressions taking 
into account measured errors on both variables. The best-fit 
relation and the residuals of the best-fit are shown in Fig- 
ure [T] The standard deviation of these residuals is 1.5 mm. 
The linear regression gives a scale factor of 0.9±0.1, with an 
offset at zero of -1.2±0.1 mm. 

The 940nm-radiometer provided the statistical values 
taken as representative of the ORM site until now. There- 
fore, we adopted the 940nm-radiometer as the reference in- 
strument to calibrate the long-term GPS time series in or- 
der to compensate the scale and offset errors introduced in 
the data processing. The standard deviation of the residu- 
als was chosen as the uncertainty associated to the PWV 
estimations from the GPS measurements at ORM. Figure 
[2] presents the full time series of PWV data derived from 
GPS measurements above ORM in the period from June 
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Figure 1. Comparison of PWV estimations from GPS and : (a)-(b) 940nm-radiometer measurements at La Palma site; (c)-(d) PWV 
estimations from atmospheric opacity at 225 GHz at Mauna Kea site, ( a) and (c) T lie solid line corresponds to the best fit derived from 
a Bayesian approach to linear regressions with errors in both variables l lKellvll2007h . (b) and (d) The residuals of the best-fit model. 



2001 to December 2008. This figure also includes the frac- 
tion of measurements for which PWV was found to be below 
a given value. 

Unlike other techniques, GPS provides PWV estima- 
tions in nearly all atmospheric conditions, including those of 
local fog, dust and/or presence of clouds (non-photometric 
conditions). These conditions will induce large PWV values 
that are not considerered in PWV time series obtained from 
940nm-radiometer measurements. 

2.2 PWV data at Mauna Kea 

Mauna Kea is an excellent astronomical site for thermal 
IR observations. The Caltech Submillimeter Observatory 
(CSO) has two monitors (Latitude 19.8224995 North, Lon- 
gitude -155.475844 West; Altitude: 4070 m above sea level) 
of atmospheric opacity, one operating at 225 GHz and an- 
other at a wavelength of 350 pim permanently installed 
since the nineties. The characterist ics of these radiome- 
ters have been extensively descri bed (|Liu Ill987l : iMcKinnon I 
Il987l : IChamberlin fc BaUvl ll994D . More information about 
these monitors is provided in the CSO Tau monitor web- 
page (http://puu oo. submm.caltech.edu/ ). Previous work 
l|Davis et al.lll997l ') has shown that the atmospheric opacity 
at 225 GHz, T22iGHz , is related to the column abundance of 
water vapor in mm, namely: 

PWV = 20 X {T225GH. - 0.016) (4) 

The time series of T225GHz measurements above Mauna 
Kea were collected from the CSO archive for the period June 
2001 to December 2008. The temporal sampling of T225GHz 
is about 10 minutes. The PWV was estimated (PWV225GHZ 
hereafter) from T225GHZ data using equation |4] Statisti- 
cal results based on this refere nce PWV225G-gz time series 
have been alr eady presented l|Garcfa- Lorenzo et al.l l2009l : 
lOtarola et aH 2010). 

A permanent GPS station (labelled MKEA) is located 
(Latitude: 19.8014 North, Longitude:-155.456 West and Al- 
titude: 3755 m above sea level) at approximately 4 km from 
Mauna Kea observatory. The coordinate estimation per pro- 



cess for the GPS antenna at MKEA produces height esti- 
mates which are very consistent day to day. This confirms 
that the station is correctly installed, stable and that the 
data is relevant in this scientific context for water vapor es- 
timations. It is important to point out that the statistical 
results from PWV derived using these GPS data might not 
be representative of the PWV conditions at Mauna Kea as- 
tronomical site as the GPS station is four km distant from 
the observatory, although we assume hereafter a similar be- 
havoiur of the PWV in both locations. The difference in alti- 
tude between Mauna Kea site and the GPS (~ 315m) might 
introduce an aditional bias in the PWV estimations derived 
from GPS relative to the values obtained from T225GHz- 
This bias could be up to 1 mm assuming a similar vari- 
ation of PWV with h eight as the found i n Atacama site s 
(iGiovanelh et al.ll200ll ) or in Pico Veleta (jOuesadal [rgsgj ). 
The data from the MKEA GPS station for the period June 
2001-December 2008 were proccessed using the same global 
GPS data procedure that we used for the GPS station at La 
Palma site. In order to compare the GPS time series esti- 
mations of PWV to PWV225GHZ , the CSO Tau monitor is 
adopted as the reference instrument. The PWV225GHZ data 
were re-sampled to the same temporal resolution of the GPS 
time series (2 hours). The standard deviation of the averaged 
PWV225Gifz in the two-hour sampling is taken as the error 
of each measurement. Errors of about 1 mm are assumed 
for the PWV derived from the proccesing of the GPS data 
(PWVgps hereafter). The number of data to be compared 
is 20880. The Pearson correlation coefficient obtained when 
comparing PWVgps and PWV225GH2 is 0.91. Several ran- 
dom selections of about 500 data were performed in order 
to apply the same Bayesian approach to linear regressions 
than that used for La Palma data. The result of the linear 
regression gives a scale factor of 0.9±0.1 and an intercept 
of -4.3±0.1 mm (Fig. [l}. This offset comes from a combi- 
nation of modeling errors in the GPS treatment and the 
difference in altitude between the GPS station and the CSO 
Tau monitor. These scale and offset factors were applied in 
order to compensate the already explained errors (section 
§2.1) introduced in the GPS data processing. The standard 
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Figure 2. (a) The two-hourly time series of PWV estimations from GPS measurements at ORM from June 2001 to December 2008. 
(b) Cumulative distribution of PWV measurements (solid line) above ORM. Dotted lines represent the uncertainties due to both the 
measurements and the calibration. 



deviation of the residuals to the linear approach is 1.6 mm, 
which it is assumed hereafter as the uncertanties associated 
to PWV derived from GPS data at Mauna Kea. The full 
time series of PWV from GPS measurements obtained for 
MKEA in the period from June 2001 to December 2008, and 
the cumulative distribution of the measurements are shown 
in Figure [3l 



3 ANALYSIS AND RESULTS 

GPS measurements provide PWV estimations almost con- 
tinuously, including daytime and nighttime data. The av- 
erage PWV at ORM was 4.9 mm, with an standard de- 
viation of 3.7 mm. The median value considering all the 
estimations in the PWV time series for ORM is 3.9 mm 
(Table [TJ. Ephemerides from the Nautical Almanac of 
the Real Observatorio de la Armada de San Fernando 
("http:/ /www.armada.mde.es ) has been used to divide the 
data in day and nighttime series. The average daytime PWV 
is 5.1±3.7 mm, whereas the average nighttime PWV is 
4.8±3.7 mm (Table [T]). The uncertainties indicate only the 
standard deviation of the averaged measurements. 

The statistical study of PWV variability at different 
temporal scales is of relevant importance for an effective 
scheduling at an astronomical observatory with IR capa- 
bilities. We have used the time series of PWV estimations 
from GPS data above ORM to study the variability of this 
parameter at seasonal, monthly, and daily scales, but it is 



important to note that weather patterns are highly variable 
on all timescales. Variability for periods larger than a year 
(e.g. associated to North Atlantic oscillation events and/or a 
global warming) can not be raised with our current 7.5-year 
time series. 



3.1 Daily variations 

It is expected that larger temperatures result in larger PWV 
(as it is expressed in equation [2|. Since temperature follows 
a clear diurnal variation (|Yagal ,2002) , it might be likely 
that PWV follows a similar trend. The two-hourly time se- 
ries above La Palma site were used to analyse the variations 
of PWV in a 24-hour period. The global statistics (Table [TJ 
suggests a slight difference between day and night values of 
PWV, with larger daytime values. However, such differences 
between daytime and nighttime seems to not follow a clear 
daily behaviour, in the sense that the minimum and maxi- 
mum PWV during a full day does not have a defined period. 
Indeed, maxima or minima values of PWV might appear at 
any hour during the day as it is illustrated in Fig. [l] The 
distribution of the maximum PWV presents two peaks, at 
around noon and midnight, while the distribution of min- 
imum PWV presents a peak only at around midnight. We 
have tested that this distribution does not change signifi- 
catly if we only take into account the days from a particular 
season or month. Therefore, a clear daily behaviour of the 
PWV is not found at ORM. 
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Time (Years) PWV (mm) 

Figure 3. (a) The two-hourly time series of PWV estimations from GPS measurements close to Mauna Kea astronomical site from June 
2001 to December 2008. (b) Cumulative distribution of PWV measurements (solid line) above Mauna Kea. Dotted lines represent the 
uncertainties due to both the measurements and the calibration. 
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Figure 4. Distribution of hours at which maximum (continuous line) and minimum (dashed line) PWV are reached during a 24-hour 
period. 
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Table 1. Global and seasonal statistical results of precipitable water vapor (in mm) above the Roque de los Muchachos Observatory. 
The percentage of time in which the PWV is within a defined range is also shown. 







Global 


Winter 


Spring 


Summer 


Autumn 


AH: 
















Mean 


4.9 


3.3 


3.5 


6.7 


5.6 




cr 


3.7 


2.3 


2.4 


4.3 


3.8 




N 


30920 


6886 


7200 


8655 


8179 




10% 


1.2 


0.9 


0.9 


2.0 


1.4 




25% 


2.3 


1.7 


1.8 


3.2 


2.6 




50% 


3.9 


2.9 


3.0 


5.7 


4.8 




75% 


6.7 


4.4 


4.7 


9.4 


7.9 




90% 


10.3 


6.3 


6.6 


13.0 


11.13 


All: 


PWV Range 














< 2 mm 


217o 


3l7o 


29% 


9%) 


17%) 




2-3 mm 


17% 


21% 


20% 


13% 


13% 




3-6 mm 


33% 


36% 


37% 


30% 


31% 




6-10 mm 


18% 


10% 


12% 


26% 


25% 




^ 10 mm 


11% 


2 % 


2% 


22% 


14% 


Day: 
















Mean 


5.1 


3.4 


3.6 


6.8 


5.7 




(7 


3.7 


2.2 


2.3 


4.2 


3.8 




N 


12967 


2514 


3223 


3904 


3326 




10% 


1.4 


1.0 


1.0 


2.2 


1.5 




25% 


2.4 


1.9 


2.0 


3.4 


2.7 




50% 


4.1 


3.0 


3.2 


5.7 


4.9 




75% 


6.8 


4.5 


4.9 


9.4 


8.1 




90% 


10.5 


6.3 


6.8 


12.7 


11.5 


Day: 


PWV Range 














< 2 mm 


18% 


28% 


26% 


8% 


16% 




2-3 mm 


17% 


22% 


20% 


12% 


14% 




3-6 mm 


35% 


38% 


39% 


32% 


30% 




6-10 mm 


19% 


10% 


13% 


26% 


25% 




> 10 mm 


11% 


2% 


2% 


22% 


15% 


Night: 
















Mean 


4.8 


3.3 


3.3 


6.7 


5.7 




a 


3.7 


2.3 


2.4 


4.4 


3.7 




N 


17921 


4340 


3977 


4751 


4853 




10% 


1.2 


0.8 


0.8 


1.9 


1.4 




25% 


2.1 


1.6 


1.6 


3.1 


2.6 




50% 


3.8 


2.8 


2.9 


5.7 


4.7 




75% 


6.6 


4.4 


4.5 


9.5 


7.8 




90% 


10.2 


6.3 


6.4 


13.1 


11.0 


Night: 


PWV Range 














< 2 mm 


23% 


33% 


32% 


11% 


17% 




2-3 mm 


16% 


20% 


20% 


13% 


13% 




3-6 mm 


32% 


35% 


35% 


29% 


31% 




6-10 mm 


18% 


10% 


11% 


25% 


25% 




> 10 mm 


11% 


2% 


2% 


22% 


14% 



3.2 Seasonal variations 

We have derived the global and seasonal statistics of PWV 
for the period June 2001-December 2008 (see Table [l]). Sea- 
sonal differences in the statistical values for daytime and 
nighttime PWV were found. The smallest statistical PWV 
corresponds to winter and spring, wheareas the largest oc- 
curs in summer. Table [1] also presents the percentage of time 
in which PWV lies on a determined range of values defined 
following the criteria t o classify the observa tional conditions 
for IR observations in iKidger et al.l l|l998l ) . Approximately 



38% of the time (including day and night) the ORM presents 
PWV values smaller than 3 mm, this percentage increases 
to more than 52% for winter and spring nighttime. On the 
contrary, values larger than 6 mm are found 29% of the total 
time, although this percentage is only about 12% in winter 
and spring. 

We have also derived the monthly average PWV for 
the different years and for the full period in order to study 
in detail the seasonal trend (figure [5]). The PWV presents 
a clear seasonal behaviour, being September and April the 
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statistically worst and best months, respectively, in terms of 
PWV. The seasonal behaviour does not significantly differ 
among the years, showing the largest PWV values during 
August-September and the lowest during March- April. The 
latter results were also found if we considered separately 
the daytime or nighttime series. This seasonal behaviour 
could be related to the seasonal variation of the temper- 
ature (see e.g. the yearly temperature graph for ORM at 
[http:/ /catserver. ing.iac.es/weather/archive/m rtg/temp_out . 
as the PWV strongly depends on this meteorological 
variable. Seasonal trends of the PWV have been also 
reported for other a stronomical sites in both hemispheres 
l|Otarola et al.ll20ld ). 

If we only take into account the nighttime series, 39% 
of the time PWV are lower than 3 mm (table [TJ . We have 
calculated the monthly percentage of night time at which 
the PWV conditions are good or excellent for IR observa- 
tions, that is PWV ^ 3 mm (Figure[6|. This percentage also 
varies along the months. In February, March and April, more 
than 60% of the nighttime presents high quality conditions 
to perform IR observations. Even in September, the worst 
month in terms of PWV, 10% of the nighttime is highly ade- 
quated to perfom IR observations. In this sense, high quality 
IR observations at ORM could be carried out at any time 
of the year, being April the best month with ~ 63% of the 
nighttime presenting PWV ^ 3 mm. 



< fFD{M) >=< fpDioo) > x[l - exp{~{At/Ty)] (6) 

The <fFD(At)> should range from <fFD(0)>=0 to a 
saturated value <fFD(oo)>. From the PWV time series for 
the ORM, we have computed the average normalized differ- 
ence for the complete data set for At ranging from to 72 
hours in steps of 2 hours following equation[5] (see Fig.[7Ka)). 
The characteristic time and the growth rate (r and 7) have 
j^^^jjbeen determined by fitting the equation[6]to <fFD(At)> us- 
ing a gradient-expansion algorithm (task CURVEFIT from 
the IDL Astronomy User's Library) and minimizing ■ The 
derived values are: r = 46.28 hours, 7 = 0.43, and < 
fFD(oo)> =0.49 mm. The PWV at a time t+At will be, 
on average, within the range: PWV(t) x [l±fFD(At)]. Let us 
assume that PWV=2.5 mm at a given time t, the PWV ex- 
pected at a later time t±At (with At ranging from to 72 
hours) will be in the range between the two curves in figure 
Elb). 

The AD method was applied bv lSkidmore et all (|2009l ) 
to study the temporal variability of the seeing at the five 
candidate sites for the Thirty Meter Telescope (TMT). The 
AD method provides information of the absolute range of 
expected PWV at a time t+At knowing the PWV at a given 
time t. The AD is calculated over the entire data set simply 
applying: 



3.3 Temporal Stability 

Temporal stabili ty is an important property in site- 
testing analysis (iRacind 1 19961 : iMufioz-Tunon et al.l 1 19971 : 
ISkidmore et al.ll2009l '). although it is not usually taken into 
account. Different astronomical sites with similar average 
values for a particular parameter could strongly differ in the 
temporal stability of that parameter. In this sense, the char- 
acterization of the temporal evolution of a parameter might 
be used to optimize the operability of astronomical facilities. 
The temporal fluctuations of PWV drastically affects the 
schedule of an observing night at a telescope/observatory 
working in queue mode. Therefore, estimations of the range 
of PWV fluctuations on a given time scale is an important 
factor to take into account for queue-scheduled IR obser- 
vations. The long-time series of PWV estimations recorded 
at ORM using a two-hour sampling of GPS measurements 
allows us to evaluate the temporal variability of the PWV. 
Variations of PWV in time scales smaller than two hours 
can not be studied with the current time series. 

In order to study the temporal variability of the PWV 
at the ORM, we have applied the fractional difference (FD) 
and the absolute difference (AD ) methods. The FD method 
was introduced bv lRaciii3 l| 19961 ) to study the temporal fluc- 
tuations of the seeing at Mauna Kea. The FD method is 
based on a function fFD(At): 



fFD{At) =< 



{PWVjt + At) - PWVit)\ 
PWVit + At) + PWV{t) 



fAD{At) =< \PWV{t + At) - PWV{t)\ > 



(7) 



>,0 ^ fpoiAt) ^ 1(5) 



which measures the characteristic range of relative 
changes as a function of the relative time At. Assuming a 
log-normal distribution of the data, the characteristic fluc- 
tuation time (tfd) and the growth rate (7F13) can be deter- 
mined by fltting the observations to the following function: 



The average <fAD {At)> over the entire data set derived 
from the PWV time series for the ORM is shown in figure 
\8la). Let us again assumed that PWV(t)=2.5 mm, the PWV 
at a time t+At should be in the range between the two 
curves in figure [SKb). 

According to figures [Tjb) and [Sjb), and considering 
At=24 hours, the PWV(t+24) should be in the range 1.86- 
3.14 mm using the FD approach or in the range 0.29-4.70 
mm when the AD method is prefered. The AD method seems 
to be more restrictive than the FD method, but the reader 
should take into account that both procedures are averaging 
over a large number of data and the dispertion is large in 
both cases. 

We have also calculated the average temporal range in 
which the PWV is always below a defined value PWVo (Fig. 
[9]). We obtained that the average temporal range in which 
conditions are good or excellent for IR observations (PWV^ 
3 mm) at the ORM is around 16.9 hours. 

The three methods that were applied to study the tem- 
poral variability of the PWV represent statistical analysis of 
such temporal variations. A steep rise of the PWV in a short 
time scale can not be ruled out. Indeed, we found that differ- 
ences larger than 1 mm in consecutive (two hours sampling) 
PWV measurements occur in ~9% of the cases. This per- 
centage decreases to ~1.5% for differences larger than 2 mm 
in consecutive measurements. A significant improvement in 
PWV might also be possible. In the two-hour PWV time 
series for ORM, improvements larger than 1 mm or 2mm 
also occurs in ~9% and ^^1.6% of the cases, respectively. 



3.4 PWV statistics at Mauna Kea derived from 
GPS data 

In order to compare the statistical results on PWV derived 
for ORM with a higher altitude site, we also carried out an 
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Figure 5. (a) Monthly statistical values of PWV for a 7.5-year period above ORM. This time series constitues a smoothed dataset of the 
two-hour time series derived from GPS measurements recorded at the ORM (including night and day data); (b) The monthly statistics 
of PWV for the period June 2001-December 2008 for ORM. In both plots, dots correspond to the average values and triangles to median 
values. Errorbars only indicate the standard deviation. 



statistical analysis of PWV estimations derived from GPS 
measurement at a location closer to Mauna Kea observatory. 
Previous statistical results of PWV h ave been alre a dy ob - 
tained from other technique s (see e.g. I Savers et al.l l|2009l ): 
iGarcfa- Lorenzo et al.l (|2009l ')') . In section tj3.5l we will sum- 
marize PWV statistics for different worldwide astronomical 
sites obtained from data recorded using different techniques. 

Table [5] presents the statistical results of the PWV time 
series derived from GPS measurements for the period from 
June 2001 to December 2008 at four km distant from Mauna 
Kea observatory. The global average PWV at Mauna Kea 
is almost 2 mm smaller than the corresponding mean PWV 
at ORM, which is not surprising due to the large difference 
in altitude between both sites (more than 1300 m). Statisti- 
cally, Summer and Autumn present slightly larger PWV at 
Mauna Kea than Winter and Spring. However, a clear sea- 
sonal trend in the PWV behaviour is not clear from figures 
[ToTa') and[TOi;b). The PWV is statistically smaller than 3 mm 
in 63% of the time (night -1-day), being this value almost in- 



dependent of the season, although during spring nighttime 
the percentage increases to 70%. PWV measurements larger 
than 6 mm appear around 12% of the time at any season. 
The average temporal range in which conditions are good or 
excellent for IR observations (PWV^ 3 mm) at Mauna Kea 
is 23.4 hours (see figure [9]). 

Comparing tables [T] and [2] we can deduce that for each 
10 hours of usefuU nighttime for IR observations (PWV ^ 
3mm) at Mauna Kea, the ORM will present 6 useful hours. 
During Summer and Autumn, only 3.9 and 4.7 hours are use- 
ful for IR observations at ORM for each 10 useful hours at 
Mauna Kea, whereas during Winter and Spring these figures 
increase to 7.9 and 7.4 useful hours, respectively, for each 10 
hours at Mauna Kea. Fair or mediocre conditions for IR ob- 
servations (3mm ^ PWV ^ 6 mm) differ in around 10% 
when comparing results from Mauna Kea and ORM, being 
in any case larger at ORM. During Winter and Spring, ORM 
presents approximately the same percentage of bad condi- 
tions (PWV ^ 10 mm) than the percentage at Mauna Kea 
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Figure 6. Monthly percentage of night time at which PWV ^ 3 mm for the period June 2001-December 2008 at the ORM. 
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Figure 7. (a) Fractional difference of the PWV measurements as a function of the time interval averaged over the full data set. The 
standard deviation of the <f^£3(A t)> is 0.22 in average, (b) PWV range expected (continua lines) at t+At when assuming a PWV(t)=2.5 
mm (long-dashed line) applying the FD method. The short-dashed line indicates the characteristic time r. 



along the seasons. Moreover, the average temporal range in 
which conditions are stably excellente or good (PWY^S 3 
mm) are similar (within uncertainties) at Mauna Kea and 
ORM. 

This comparison demostrates that excellent conditions 
for IR astronomical observations, in terms of percentage 
and stability, are also possible in sites at relative low al- 
titude (over 2000 m above sea level). This result might be 



related to the role of the troposphere thickness instead of 
the site a ltitude in the IR qual i ty of a particular astronom- 
ical site l|Garcfa- Lorenzo et al.l [20041 ). Indeed, the seasonal 
behaviour of the PWV at ORM follows a quite similar be- 
havi our than the troposphere thic kness at this site (see figure 
5 in lGarci'a-Lorenzo et al.l (|2004l )). 
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Figure 8. (a) Absolute diflference of the PWV measurements as a function of the time interval averaged over the full data set. The 
average standard deviation of the data is 2.30 mm. (b) Expected range of PWV at any time t+At assuming a PWV equal to 2.5 mm 
(dashed line) at time t using the AD approach. 
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Figure 9. (a) Average temporal interval in which the PWV is lower than a defined PWVq value at ORM; (b) Average temporal interval 
in which the PWV is lower than a defined PWVo value at Mauna Kea. 



3.5 Summary of PWV statistical results at 
astronomical sites 

The comparison of PWV statistical results between astro- 
nomical sites is very complex due to the variety of techniques 
and procedure normally used. Different techniques provide 
different temporal coverages and samplings. Moreover, some 
of the techniques used to evaluate the water vapor content 
are strongly affected by rain and clouds (e.g. radiosondes 
and WVR). Others are only applicable near the zenith (e.g. 
radiometers), whereas GPS can provide continuous PWV 
estimates depite the effects of rain, dust or clouds. In spite 
of such differences, we summarize in table |3]statistical PWV 
results for different astronomical sites located at different al- 
titudes above sea level. In the case of Mauna Kea, we have 
found median PWV values derived from four different tech- 



niques, showing significant discrepancies between the values. 
Even using the same database but different analysis period, 
statistical values may differ significantly. These differences 
also occur in the case of ORM, where Infrared sky radiance 
gives a significant lower PWV median than the other tech- 
niques used there. In both sites, GPS provides the largest 
PWV median value, as it is expected from the fact that GPS 
measures continuosly and under all weather conditions. 

In order to calculate an approximate PWV median 
value comparable to that derived from the techniques 
affected by clouds, we have selected the GPS data recorded 
during photometric nights (clear nights from clouds and/or 
aerosol affecting astronomical observations). The atmo- 
spheric extinction coefficient - a measurement of the sky 
transparency - is measured at the ORM in the V (551 nm, 
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Table 2. Global and seasonal statistical results of precipitable water vapor (in mm) above Mauna Kea region including the percentage 
of time presenting a PWV value in a determined range. 







Global 


Winter 


Spring 


Summer 


Autumn 


AH: 
















Mean 


3.0 


2.8 


2.7 


3.1 


3.2 




cr 


2.4 


2.50 


2.2 


2.2 


2.6 




N 


28607 


5541 


6897 


8417 


7752 




10% 


0.6 


0.3 


0.5 


0.8 


0.6 




25% 


1.2 


0.9 


1.2 


1.5 


1.3 




50% 


2.3 


1.9 


2.1 


2.6 


2.3 




75% 


4.1 


4.1 


3.7 


4.2 


4.3 




90% 


6.5 


6.7 


5.9 


6.2 


7.4 


All: 


PWV Range 














< 2 mm 


44% 


52%i 


4770 


36% 


43% 




2-3 mm 


19% 


13% 


19% 


21% 


19% 




3-6 mm 


25% 


22% 


24% 


32% 


23% 




6-10 mm 


11% 


12% 


9% 


10% 


13% 




^ 10 mm 


1% 


1% 


1% 


1% 


2% 


Day: 
















Mean 


3.1 


2.8 


2.8 


3.3 


3.3 




(7 


2.4 


2.5 


2.2 


2.2 


2.6 




N 


16562 


3072 


4418 


5141 


3931 




10% 


0.7 


0.4 


0.7 


0.9 


0.7 




25% 


1.3 


0.9 


1.3 


1.6 


1.3 




50% 


2.4 


2.0 


2.2 


2.7 


2.4 




75% 


4.2 


4.1 


3.8 


4.3 


4.5 




90% 


6.6 


6.7 


5.9 


6.3 


7.5 


Day: 


PWV Range 














< 2 mm 


41% 


50% 


44% 


33% 


41% 




2-3 mm 


19% 


14% 


21% 


22% 


19% 




3-6 mm 


27% 


23% 


25% 


34% 


24% 




6-10 mm 


11% 


12% 


9% 


10% 


14% 




> 10 mm 


2% 


1% 


1% 


1% 


2% 


Night: 
















Mean 


2.9 


2.7 


2.5 


2.9 


3.1 




a 


2.4 


2.5 


2.3 


2.2 


2.6 




N 


12043 


2462 


2484 


3276 


3821 




10% 


0.5 


0.3 


0.4 


0.7 


0.6 




25% 


1.1 


0.8 


0.9 


1.4 


1.2 




50% 


2.1 


1.8 


1.8 


2.4 


2.2 




75% 


3.9 


3.9 


3.5 


4.0 


4.1 




90% 


6.5 


6.7 


5.7 


5.9 


7.3 


Night 


PWV Range 














< 2 mm 


48% 


55% 


53% 


42% 


45% 




2-3 mm 


17% 


12% 


17% 


19% 


19% 




3-6 mm 


23% 


20% 


21% 


29% 


22% 




6-10 mm 


10% 


12% 


8% 


9% 


12% 




> 10 mm 


2% 


1% 


1% 


1% 


2% 



Kv hereafter) and r' (625 nm) bands since 1984 by the 
Carlsberg Automatic Meridian Circle (CAMC). We have 
downloaded the atmospheric extinction coefficient for the 
period June 2001-December 2008 from the CAMC archive 
I http://www.ast.cam.ac.uk/dwe/SRF/camc_extinction.htmll. 
We selected as photometric nights those satisfying the fol- 
lowing requirements: (i) the total number of observing hours 
is larger than seven; (ii) the total number of photometric 
data taken is larger than 90% of the total number of 
observing hours; and (iii) the average along the night 



is smaller than 0.2. 55% of the nights from June 2001 to 
December 2008 in the CAMC database follow the selection 
criteria. PWV measurements derived from GPS data during 
photometric nights were then used to derived the statistical 
values of PWV (see table [H). 

Statistical values (average and median) of PWV 
strongly improve when specific weather conditions are se- 
lected, which is usually imposed by the techniques used to 
retrive PWV values. Therefore, the comparison of PWV val- 
ues from different sites is only possible when the same tech- 
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Figure 10. (a) Monthly average PWV values for a 7.5-year period above Mauna Kea. This time series constitues a smoothed dataset 
of the two- hour time series derived from GPS data close to Mauna Kea Observatory (including night and day data); (b) The monthly 
statistics of PWV for the period June 2001-December 2008 for Mauna Kea. In both plots, dots correspond to the average values and 
triangles to median values. Errorbars only indicate the standard deviation of the data; (c) Monthly percentage of night time at which 
PWV ^ 3 mm at Mauna Kea for the period from June 2001 to December 2008. 



Table 3. Statistical results of precipitable water vapor for different astronomical sites worldwide derived from different techniques and 
procedures. 



Site 


Location: 


Height 


Median 


Technique 


Temporal 


Reference 




Lat, Lon 


(m) 


PWV (mm) 




Range 




Las Campanas 


29.01° 


S, 42.18° W 


2200 


2.8 


225GHz-radiometer 


2005/07-08 


Thomas-OsiD et al. C2007) 


Cerro Tolar 


21.96° 


S, 70.10° W 


2290 


4.02 


GOES-8 satellite 


1993/06-1996/02 


Otarola et al. ("2010) 










4.7 


Surface PWV data 


01/2004-12/2007 


Otarola et al. f2010~) 


ORM 


28.77° 


N, 17.88° W 


2395 


3.9 


940nm- Radiometer 


1996-1998 


Kideer et al. (19981 










3.9 


GPS 


2001/06-2008/12 


this work 










2.6 


IR Sky Radiance 


2000-2002 


Pinilla ("2003) 


La Silla 


29.25° 


S, 70.73° W 


2400 


3.9 


IR Sky Radiance 


1983-1989 


Sarazin fl990~) 


Paranal 


24.63° 


S, 70.40° W 


2635 


2.3 


IR Sky Radiance 


1983-1989 


Sarazin (1990) 


San Pedro Martir 


31.04° 


N, 115.47° W 


2830 


2.63 


GOES-8 satellite 


1993/06-1996/02 


Otarola et al. (2010) 










3.4 


210GHz-radiometer 


2006 


Otarola et al. (2010) 


Pico Veleta 


37.07° 


N, 3.37° W 


2850 


2.9 


940nm- Radiometer 


1984-1987 


Quesada (1989) 


Cerro Armazones 


24.58° 


S, 70.18° W 


3064 


2.87 


GOES-8 satellite 


1993/06-1996/02 


Otarola et al. (2010) 










3.2 


Surface PWV data 


01/2004-12/2007 


Otarola et al. (2010j 


Dome C 


75.06° 


S, 123.23° E 


3233 


0.34 


Satellite &; Model 


2008 


Saunders et al. (2009) 


Mauna Kea 


19.83° 


N, 155.47° W 


4205 


1.7 


225GHz-radiometer 


2001/06-2008/12 


Garcfa-Lorenzo et al. (2009) 










1.2 


Radiosondes 


1983 


Belv (1987) 










1.86 


GOES-8 satellite 


1993/06-1996/02 


Otarola et al. (2010) 










2.1 


225GHz-radiometer 


01/2004-12/2007 


Otarola et al. (2010) 










2.3 


GPS 


2001/06-2008/12 


this work 


Ridges A 


81.5° S, 73.5° E 


4053 


0.21 


Satellite k. Model 


2008 


Saunders et al. (2009) 


Dome A 


80.73° 


S, 77.3° E 


4083 


0.23 


Satellite &; Model 


2008 


Saunders et al. (2009) 


Cerro Tolonchar 


23.93° 


S, 67.97° W 


4480 


1.7 


GOES-8 satellite 


1993/06-1996/02 


Otarola et al. (2010) 










1.8 


Surface PWV data 


01/2004-12/2007 


Otarola et al. (2010) 


Chajnantor 


23.02° 


S, 67.45° W 


5080 


1.0 


Radiosondes 


1998/10 & 2000/08 


Giovanelli et al. (2001) 










1.2 


225GHz-radiometer 


1995/04-2000/04 


Thomas-OsiD et al. (2007) 
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Table 4. Seasonal statistical results of precipitable water vapor (in mm) above the Roque de los Muchachos Observatory derived from 
PWV estimations during photometric nights. 



Global Winter Spring Summer Autumn 

Night: 



Mean 


3.9 


3.2 


3.1 


4.7 


4.5 


a 


2.6 


1.9 


1.9 


2.9 


2.8 


N 


6146 


1256 


1517 


2081 


1292 


10% 


1.2 


0.9 


0.9 


1.7 


1.4 


25% 


2.1 


1.7 


1.7 


2.6 


2.3 


50% 


3.4 


2.8 


2.9 


3.9 


4.0 


75% 


5.2 


4.2 


4.2 


6.2 


6.2 


90% 


7.3 


5.7 


5.7 


8.4 


8.3 



nique and procedure have been used to adquire and analyse 
the data. 



3.6 Infrared bands and PWV 

ORM has been an astronomical site traditionally dedi- 
cated to optical and near-IR observations (from 0.4 to 2.5 
^m). The 10 meters Gran Telescopio de Canarias will ex- 
tend this range to the mid-IR, with CanariCam, a cam- 
era and spectrograph working in the thermal infrared be- 
tween ~7.5 and 25 pm. In this section, we explore the im- 
plications in terms of atmospheric transmission of PWV 
values in the standard near and mid-IR windows in as- 
tronomy. We have modelled the theoretical near and mid- 
IR transmission spectrum for La Palma site using the 
ATRAN modelling software (Lord, S.D. 1992, NASA Tech- 
nical Memor. 103957) throughout the web-based form at 



cover, etc) in the quality of dat a recorded in mid-IR bands 
was illustrated in iMason et al. using real data ob- 



http://atran.sofia.usra.edu/cgi-bin/atran/atran.cgi . We se- 



lected the closest latitude to the site allowed (30 de- 
grees), the altitude of the site (2400 m), scanning the 
PWV from 1 to 10 mm. The transmission spectrum de- 
rived from each PWV value was integrated in the spec- 
tral range covering the standard filters windows and these 
values were refered to the integrated transmission asum- 
ming a PWV equal to 1 mm (see table [5]). The the- 
oretical transmission spectrum for Mauna Kea was ex- 
tracted from a model available at the GEMINI Observatory 
web pages (http://staff.gemim.edu/~kvolk/linkpage.html) 
which is also based on the ATRAN modeling software, al- 
though it is limited in PWV values. We followed the some 
procedure than in the case of La Palma, integrating the 
transmission spectrum in each filter window and refering 
the value to the derived for 1mm of PWV. 

In both cases, similar results were found. The atmo- 
spheric transmission in bands H, K, and N decrease less 
than 5% for PWV variations from 1 to 10 mm. The atmo- 
spheric transmission in windows J, L, and M slowly decrease 
when PWV increase (up to 21% from 1 to 10 mm). The most 
significant effect occurs in band Q, where atmospheric trans- 
mission is around 50% smaller when PWV increase from 1 to 
5 mm. When PWV is 10 mm, the atmospheric transmission 
is almost 35% of the transmission with 1 mm of PWV in Q- 
band. It is also important to note that high levels of water 
vapor also increase the thermal IR background, which is one 
of the major factors limiting the atmospheric transparency 
in mid-IR observations from ground-based sites. The effects 
of observing conditions (water vapor column, airmass, cloud 



tained in Mauna Kea. 



4 SUMMARY AND CONCLUSIONS 

We have analised the water vapor content for the period 
from June 2001 to December 2008 above the ORM using 
PWV estimations from GPS data. We have verified the con- 
sistency of 940nm-radiometer and GPS estimation of PWV, 
removing the offset between both techniques. We have also 
presented statitical results for Mauna Kea for the same 
period, analysing the GPS close to Mauna Kea site with 
T225GHZ measurements. Our main results and conclusions 
may be summarized as follows: 

(i) The nighttime median PWV above ORM is 3.79 mm, 
with slightly differences between day and nighttime statis- 
tics. 

(ii) The PWV presents a clear seasonal variation at the 
ORM. Winter and spring nights present the lower PWV 
statistics, with a median value of 2.82 mm and 2.86 mm, 
respectively. 

(iii) More than 60% of the nighttime during February, 
March and April present good or excellent conditions in 
terms of water vapor (PWV ^ 3 mm) at the ORM. 

(iv) Comparing PWV statistical results for ORM and 
Mauna Kea, we deduce that for 10 hours of good condi- 
tions for IR observations at Mauna Kea during winter, the 
ORM will present 7.9 hours showing excellent conditions. 

(v) The average temporal range presenting good or excel- 
lent conditions for IR observations (PWV ^ 3mm) is compa- 
rable at ORM and Mauna Kea (16.9 hours and 23.4 hours, 
respectively). 

(vi) The comparison of PWV at different astronomical sites 
is difficult because there is not an unique defined tech- 
nique/procedure to estimate the PWV. 

GPS is a promising technique to unify the PWV esti- 
mations at many astronomical sites. 
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Table 5. Atmospheric transmission as a function of PWV relative to the atmospheric transmission extracted for 1 mm of PWV for the 
standard near and mid infrared bands. 



Near 


PWV 


Filter J 


Filter H 


Filter K 


Filter L 


Filter M 


Filter N 


Filter Q 


La Palma 


(mm) 


(1.1-1.4 iim) 


(1.5-1.8 ^lm) 


(2.0-2.4 /^m) 


(3.0-4.0 /im) 


(4.5-5.1 tim) 


(10-12 /im) 


(17.5-22.5 /jm) 




1.5 


97.4 


99.7 


99.8 


97.5 


97.5 


99.8 


89.3 




3.0 


92.5 


98.8 


99.1 


92.2 


91.8 


99.4 


68.2 




5.0 


88.7 


98.1 


98.4 


87.9 


87.2 


99.0 


54.7 




6.5 


86.8 


97.6 


97.9 


85.6 


84.9 


98.8 


48.9 




7.6 


85.8 


97.4 


97.7 


84.4 


83.9 


98.7 


47.2 




10.0 


83.2 


96.7 


96.9 


80.9 


79.4 


98.3 


36.5 


Mauna 


PWV 


Filter J 


Filter H 


Filter K 


Filter L 


Filter M 


Filter N 


Filter Q 


Kea 


(mm) 


(1.1-1.4 fim) 


(1.5-1.8 fim) 


(2.0-2.4 ^m) 


(3.0-4.0 fira) 


(4.5-5.1 ^J.m) 


(10-12 ^J.m) 


(17.5-22.5 /tm) 




1.5 


98.6 


99.8 


99.9 


99.0 


98.8 


99.1 


84.8 




3.0 


91.1 


98.7 


99.1 


93.3 


92.2 


98.2 


68.3 




5.0 


87.2 


97.9 


98.4 


89.7 


88.1 


96.9 


55.9 




6.5 


86.1 


97.6 


98.6 


89.0 


88.0 







Mauna Kea on the island of Hawaii. Both GPS stations are 
part of the international network of GPS stations EUREF 
(www.euref.eu). 

Measurements of the 225 GHz optical depths 
at Mauna Kea were recorded from the web page 
of the Caltech Submillimeter Observatory (CSO) 
I http://puuoo.submm.caltech.edu/ ). The coefficient- 
extinction data were downloaded from the database 
(http:/ /www.ast. cam.ac.uk/~dwe/SRF/camc_extinction. html 
of the Carlsberg Meridian Circle of the Isaac Newton Group 
on. 
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